Surface charge of a flat superconducting slab in the Meissner state 
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The electrostatic potential in the flat superconducting slab is evaluated in the framework of the 
Ginzburg-Landau theory extended by Bardeen to low temperatures. For magnetic fields below B c \, 
we discuss the formation of a surface charge induced by the Bernoulli potential of the suppercurrents. 



I. INTRODUCTION 



II. EXTENDED GINZBURG-LANDAU THEORY 



A diamagnetic current in a superconductor creates an 
electrostatic potential known as the Bernoulli potential 
Early studies based either on the electrodynamics 
of the ideal charged fluid jjj , on the London approxima- 
tion Q or on the local approximation of the BCS theory 
|^-|| provide the electrostatic potential as a function of 
the condensate velocity. This approach is not applicable 
to strongly inhomogoneous systems like Abrikosov vor- 
tices or surfaces where non-local quantum corrections are 
inevitable. 

A simple non-local theory of the electrostatic poten- 
tial in stationary superconductors has been introduced in 
H . It is based on the Ginzburg-Landau (GL) theory and 
neglects the so called quasiparticle screening of van Vi- 
jfeijken and Staas and the thermodynamic correction 
predicted by Adkins and Waldram || and evaluated by 
Rickayzen Q]. 

A theory of Ginzburg-Landau type which covers the 
quasiparticle screening and the thermodynamic correc- 
tions, has been introduced in Its implementation to 
the Abrikosov vortex lattice is in preparation. In this 
paper we employ a non-local approach to discuss the for- 
mation of the surface charge. We study the simple system 
of a superconducting flat slab in a parallel magnetic field 
below B c \. 

The formation of a surface charge induced by the 
Bernoulli potential has been already studied by Jake- 
man and Pike |J. Starting from the time-dependent 
Ginzburg-Landau theory, they have derived a Poisson 
equation for the electrostatic potential with screening 
over the Thomas-Fermi length and with a source term 
that includes the quasiparticle screening of van Vijfei- 
jken and Staas Except for the Thomas-Fermi screen- 
ing their approach is local, what implies that a surface 
charge is formed on the scale of the Thomas- Fermi screen- 
ing length. Here we show that the non-local correc- 
tions to the electrostatic potential are dominant over the 
Thomas- Fermi screening so that the surface charge forms 
on the scale of the GL coherence length. 



Our treatment of the electrostatic potential starts from 
the free energy 
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The kinetic energy of the condensate is of the GL form 
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where ij) is the GL wave function, m* = 2m and e* = 2e 
are mass and charge of a cooperon. The condensation 
energy follows the two-fluid model of Gorter and Casimir 
& 
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with the condensation energy per volume given by Som- 
merfeld's 7 as e con = rjT 2 . The total density is the sum 
of the condensate and normal electrons, n = 2\i/j\ 2 + n n . 
The internal energy, U, the linear coefficient of the spe- 
cific heat per volume, 7, and the critical temperature, T c , 
are functions of the density n. 
The magnetic free energy, 
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depends on the applied magnetic field B a and the internal 
field B = V x A. The Coulomb energy is treated in the 
non-relativistic approximation, 
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where p = e*\ip\ 2 + en n + pxatt is the charge deviation 
from neutrality. We do not assume any external electric 
field. 

Except for the Coulomb potential, this free energy has 
been introduced by Bardeen |l3] soon after the paper of 
Ginzburg and Landau ]lO| as its simple extension to low 
temperatures. 
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In the stationary state the free energy has a minimum 
which can be found by the variation principle Jl3[ . Vari- 
ations with respect to tp, A and n„ lead to a set of ex- 
tended Ginzburg-Landau equations: 
the Schrodinger equation, 
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the Maxwell equation, 
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m* 

and the screened Poisson equation, 
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We note that variations result directly in a condition for 
the charge density. The electrostatic potential tp has been 
introduced to simplify expressions. It links to the charge 
density via the plain Poisson equation 
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The last term on the right hand side of (||) is the non- 
local Bernoulli potential ||. It is proportional to \tp\ 2 /n 
in the spirit of the quasiparticle screening |Q . The two 
other terms combine into the thermodynamic correction 
introduced by Rickayzen Q . Note that the second term is 
the thermoelectric field of the normal metal reduced by 
a factor y/l — 2\ip\ 2 /n. The screening on the Thomas- 
Fermi length is covered by the differential term in the 
left hand side of (||). 

The usual GL theory has two independent variables, 
A, tp and consequently only two differential equations 
have to be solved. Including the electrostatic potential 
we have four indeterminates A, tp, tp, n n so that the set 
of four coupled equations (|| - |9|) should be solved. In 
reality the problem is not so complicated. The second 
GL equation ffl) is independent of n n and tp. The first 
GL equation (fm is independent of tp while it depends on 
n„ exclusively via n — n n + 2\tp\ 2 . Since n is very large 
in metals and charge perturbations are rather small, one 
can neglect deviations from the charge neutrality in the 
first GL equation. In this approximation, equations (|6]|7]) 
are solved separately. Once A and ip are known, equa- 
tion (||) provides us with the electrostatic potential and 
the charge is found from the Poisson equation (p|). 



III. FLAT SLAB IN MAGNETIC FIELD 

Let us demonstrate the behavior of the electrostatic 
potential in a superconducting slab in parallel magnetic 
field. We assume a superconductor of type II in a mag- 
netic field below B c \, i.e., in the Meissner state with the 
magnetic field penetrating from the surface. 

We take the direction x perpendicular to the slab lim- 
ited to the interval (— d, d). The magnetic field we choose 
in the direction z, B = (0, 0, B), and the vector potential 
points in the direction y, A = (0, A, 0). All functions de- 
pend exclusively on the coordinate x and the wave func- 
tion is real. 

We rescale all functions into their dimensionless coun- 
terparts, x = x\o, A = Jl$o/(27tAo), tp — tpy^n/2, 
etp = tph /(AmX 2 )), where the London penetration depth 
at zero temperature is given by Aq = m/(e 2 nno). Note 
that our scale does not change with the temperature. In 
the expressions below we skip tildes denoting new func- 
tions. Equations (H,0) now read 
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with S = k 2 (1 - i 2 ) 2 (l + t 2 ). As usual in the GL theory, 
the only material parameter relevant after rescaling is the 
GL parameter k defined at T c . In the Bardecn extension 
the rescaled equations also depend on the reduced tem- 
perature t — T/T c . 

As the boundary condition we use that the vector po- 
tential is anti-symmetric, A(—x) = —A{x), with the 
value of the derivative at the surface given by the ap- 
plied magnetic field B a . The wave function has a zero 
derivative at the surface. 

The third GL equation (ra) in the dimensionless units, 
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depends on two material parameters, 
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For Niobium these parameters can be derived from the 
McMillan formula jlij and chemical trends Jlj| giving 
values Ci = 1.9 and C 2 = 0.42. 

The reduced Thomas-Fermi screening length, 
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is very small. For Niobium Atf = 0.7 A and Ao = 390 A. 
The Thomas-Fermi screening thus enter equation ( p^ ) 
with factor r 2 = 3 • 10~ 6 . With accuracy of the order 
of 1CP 5 , the solution of ( |l2| ) is ip = ip- m + <pf Tee , with the 
induced term, 
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and the free solution 
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The free solution decreases near the surface on the scal< 
of the Thomas-Fermi screening length r = 1.8 • 10 -3 . 

The amplitude of the free solution has to be selected sc 
that the slab remains charge neutral. In other words, th< 
free solution supplies the surface charge which is localizec 
on the Thomas-Fermi screening length. The neutrality o 
the slab is equivalent to the condition that the electric 
field vanishes at the surface. The amplitude of the fret 
solution is thus given by the condition dtp/dx = on th< 
surface. 

Following it is possible to show that the free solu 
tion is zero. The induced potential ip in has zero derivativi 
on the surface as it follows from (|l6|) and the GL bound 
ary condition, dtp/dx = 0. Accordingly, in the non-loca 
approach one has tpo = 0. 



IV. LOCAL VERSUS NON-LOCAL PICTURE 

The local approximation corresponds to a neglect o 
the gradients in the GL equation (11). The density o 
the condensate, ip 2 , is thus found from a simple loca 
condition, 
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From ( |l8|) one obtains ip 2 as a function of the vector po 
tential A. This function is first used in (jlTj) to solve for / 
and finally substituted into ( |l2|) to find the electrostatic 
potential ip. 

The wave function ip found within the local approxi 
mation ([TJ|) does not satisfy the GL boundary conditioi 
dtp/dx — 0. Accordingly, the derivative of the inducec 
potential is non-zero, ip- m ^ 0, so that the free solutioi 
has to be supplemented to maintain the charge neutral 
ity. 

Now we compare the full solution of (10 - |l^) wit! 
the local approximation. In our treatment we assumt 



Niobium with Oxygen impurities of concentration giving 
k = 1.5. The width of the slab we take 6Ao so that the 
magnetic field is screened from the bulk at low tempera- 
tures but penetrates the whole slab close to T c . 

At low temperatures one finds that the magnetic field 
suppresses the condensate only in a negligible fraction. 
Since the condensate density is nearly constant, it follows 
from ( fTo| ) that the magnetic field differs only negligibly 
from the London approximation, B = Bq cosh(a;/A), with 
Bo = B a j cosh(d/A) and Awl. Of course, since cor- 
rections beyond the London approximation are negligi- 
ble, the local and non-local approaches yield very similar 
magnetic fields. 




FIG. 1. The electrostatic potential at t — 0.1 and B w B c i. 
The local approximation (dashed line) differs from the full 
non-local solution (full line) mainly close to the surface shown 
in the insert. 

In figure 1 we show the electrostatic potential. The 
local approximation agrees with the non-local approach 
very well except for the close vicinity of the surface. This 
region is affected by the GL boundary condition. The 
presented result is for a magnetic field close to B c \ . 
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FIG. 2. The charge distribution corresponding to the elec- 
trostatic potential from Fig. 1. 

The formation of the surface charge is shown in Fig. 2. 
The depleted charge appears in the surface layer of char- 
acteristic width which is small compared to the London 
penetration depth or the GL coherence length, however 
still very large on the scale of the Thomas-Fermi screen- 
ing. 
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For temperatures close to T c the local approximation 
becomes unreliable if one looks for the electrostatic po- 
tential or the related charge distribution. This contrasts 
with the magnetic properties. For temperature t = 0.9 
numerical results confirm that for a weak applied field, 
B = 0.1 B c i, the magnetic field is quite well described 
by the London approximation, B = Bq cosh(x/A), with 
Bo = B a j cosh(d/A) and A = 1.7. It again follows from 
a small effect of the magnetic field on the wave function 
which is nearly constant across the slab and reduced by 
5% compared to the bulk value ipo — \J\ — t^ = 0.6, see 
Fig. 3. One can see that the local approximation for the 
deviation of the wave function is rather bad in this case. 
Indeed, the coherence length is comparable to the width 
of the slab so that the GL boundary condition is essential 
in the whole volume. 




FIG. 4. The electrostatic potential (p = ipi n + iflf iee for 
t — 0.9 and B = 0.1B c i. The constant term is chosen so 
that both potentials (16) and (17) would reach zero in the 
bulk of a thick sample. 



FIG. 3. Deviation of the wave function rp from its equilib- 
rium value, ipo — \fl — t 4 w 0.6 for t — 0.9. Already for a 
weak field, B = 0.1 B c \, the local approximation (dashed line) 
is very different from the full non-local solution (full line). 
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FIG. 5. The charge distribution corresponding to the elec- 
trostatic potential from Fig. 4. 

The local and non-local values of the electrostatic po- 
tential shown in Fig. 4 differ appreciably. First, the non- 
local potential is shifted down compared to the local ap- 
proximation. We note that in both approximations the 
initial is identified with the bulk of a thick sample. The 
constant shift, however, is not important for the distri- 
bution of the charge shown in Fig. 5. One can see that 
in the local approximation the charge is positive every- 
where. The charge neutrality is maintained by an invis- 
ibly thin depleted layer on the scale r w 10~ 3 . In the 
non-local approach, the layer of the depleted charge ex- 
tends over the region comparable with the GL coherence 
length. 

In conclusion, the formation of the surface charge due 
to the electrostatic potential caused by diamagnetic cur- 
rents has been treated within the local and non-local ap- 
proaches. While the local approach requires to include 
the surface on the scale of the Thomas-Fermi screening 
length, the non-local approach predicts that the surface 
charge extends over a layer of width comparable with the 
coherence length. 
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